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1.  Introduction  and  Background 

Statistical  data  from  household  surveys  and  administrative 
records  are  widely  used  in  formulas  to  allocate  grant-in-aid  funds  or 
other  governmental  benefits  to  state  and  local  jurisdictions.  A  review 
of  U.S.  programs  appears  in  Gonzalez  (1978).  An  explicit  tradeoff  be¬ 
tween  data  collection  costs  in  repeated  surveys  and  an  imputed  cost  of 
misallocation  of  resources  due  to  the  use  of  imprecise  survey  data  in 
setting  allocation  levels  was  considered  in  Smith  and  Zalkind  (1978)  as 
a  scalar  problem  in  deterministic  inventory  theory.  The  present  paper 
places  such  problems  in  a  control  theory  framework  accommodating  vector 
linear  models  of  nonstationary  economic  processes  measured  by  noisy 
multi-item  repeated  surveys.  Optimal  rules  are  found  for  control  of  a 
sequence  of  surveys,  where  there  is  a  fixed  charge  plus  unit  costs  of 


*An  earlier  version  of  this  paper  was  presented  at  the  2nd 
Economics  and  Control  Conference,  Princeton  University,  Princeton, 
New  Jersey,  June  2-4,  1980. 


T-428 


surveying  and  an  imputed  quadratic  loss  associated  with  the  imprecision 

! 

of  the  resulting  state  estimates  of  the  observed  economic  process. 

In  this  paper  we  analyze  one  model  in  which  the  number  of  com¬ 
pleted  interviews  in  a  survey  is  equal  to  the  designated  number  of 
interviews  and  two  models  in  which  the  number  of  completed  interviews  is 
a  random  variable.  The  models  are  analyzed  with  the  aid  of  an  equiva¬ 
lent  sample  size  form  of  the  Kalman  filter  that  was  derived  in  Smith 
(1979).  A  related  inventory  analysis  of  single-item  repeated  surveys 
under  random-yield  sampling  is  set  forth  in  Smith  and  Zalkind  (1980). 

Methods  for  time  series  analysis  of  repeated  surveys  are  presented  in 
Sct.Lt  and  Smith  (1974).  Problems  of  the  choice  of  loss  function  for 
misallocation  error  are  of  practical  importance,  but  they  are  not 
treated  in  the  present  paper,  nor  are  problems  of  "distributive  equity" 
in  allocations  across  states  or  local  jurisdictions,  as  discussed  in 
Spencer  (1979).  In  the  models  of  this  paper,  we  are  estimating  (track¬ 
ing  but  not  controlling)  the  state  of  the  observed  economic  process, 
and  we  are  controlling  the  survey  measurement  subsystem;  see  Meier, 

Peschon,  and  Dressier  (1967).  Also  see  Aoki  and  Li  (1968),  which  dis¬ 
cusses  control  problems  with  cost  for  observation,  although  not  treating 
the  fixed-plus-variable  cost  structure  of  the  present  paper. 

2 .  Formulation  of  the  Models 

In  this  section  we  formulate  the  models  to  be  discussed  in  the 
present  paper.  In  the  first  two  models  (Sections  3  and  4)  we  are  inter¬ 
ested  in  the  evolution  of  x(j)  ,  a  multivariate  socioeconomic  process 
(e.g.,  money  income  or  proportion  of  a  population  in  poverty),  through 
discrete  time  points  j  ,  j=0,l,2,...  .  Here  x(j)  is  an  mxl  state 
vector,  and  we  assume  that  it  evolves  according  to  the  vector  random 
walk 

x(j+l)  *  x(j)  +  w(j+l)  ,  (1) 

where  w(j)  -  N(0,Q)  .  That  is,  w(j)  has  a  multivariate  normal  dis¬ 
tribution  with  a  zero  mean  vector  and  known  process  disturbance 
covariance  matrix  Q  .  Surveys  of  the  process  x(j)  are  taken  once 
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every  T  units  of  time  at  k  =  T  ,  2T,  3T,  rT,  ...  .  (In  one  of 

our  stochastic  models,  T  may  vary.)  We  also  assume  that  x(0)  ~ 
N(E[x(0)]  ,  C(0))  ,  where  E[x(0)]  and  C(0)  are  known. 

At  each  survey  time  n^(k)  individuals  are  sampled,  but  only 

N(k)  completed  interviews  are  obtained.  We  assume  that  P[N(k)>l]  =  1 
The  difference  n^(k)  ~  N(k)  is  due  to  refusals  to  respond  or  incom¬ 
plete  observations.  We  assume  that  y(k)  ,  the  m-dimensional  survey 
measurement,  is  represented  by  the  survey  equation 

N(k) 

y(k)  =  x(k)  +  (l/N(k))  l  u  (k)  ,  (2) 

i*0  ~ 

where  each  u.,(k)  is  an  mxl  vector  that  is  distributed  normally 

with  a  zero  mean  and  measurement  noise  covariance  matrix  R  .  The  ma¬ 
trix  R  is  assumed  to  be  known,  time  invariant,  symmetric,  and  positive 
definite.  Both  w(j)  and  the  u^(k)  are  assumed  to  be  serially  uncor¬ 
related. 

We  may  find  the  best  linear  estimate  of  x(rT+j)  ,  given 

E[x(0)]  ,  C(0)  ,  and  the  r  measurements  yr  =  (y(T),  y(2T),  ..., 

r 

y(rT))  ,  as  a  linear  combination  of  the  observations  y  .  We  define 

0  1 

the  estimation  error  x  (k+j [ k)  as 

x£(k+j|k)  =  x(k+j )  -  x(k+j | k)  , 

where 

x(k+j |k)  =  E [ x ( k+j ) jyr]  . 

We  also  define  the  j-step-ahead  estimation  error  covariance  matrix 
C(k+j | k)  as  the  conditional  expectation  of  the  outer  product  of  the 
j-step-ahead  estimation  error  vector  given  the  r  measurements;  namely, 

C(k+j|k)  =E[xe(k+j|k)  (xe(k+j | k) ) ’ ]  .  (3) 


-  3  - 
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Wc  now  state  a  general  optimal  filter  theorem  for  the  vector 
random  walk  models.  Related  proofs  and  derivations  of  Kalman-type  fil¬ 
ters  for  vector  models  may  be  found  in  Jazwinski  (1970),  Melsa  and 
Cohn  (1978),  or  Sage  and  White  (1977). 


Optimal  Filter  Theorem  for  Vector  Random  Walk  Models .  The  optimal  (min¬ 
imum  variance)  filter  for  the  discrete  system  (i),  (2)  consists  of  dif¬ 
ference  equations  for  the  conditional  mean  x(k+jjk)  and  the  estimation 

error  covariance  matrix  C(k+i!k)  ,  for  T=1,2,...,T  : 

'  max 


he  tween  surveys , 

x(k+j|k)  =  x(k]k)  ,  for  k=T,2T,...  ,  and  j=l,...,T  , 
C(k+j | k)  =  C(k| k)  +  jQ  ; 


(4) 


at  surveys 

x(k | k)  =  x(k  j  k-T)  +  K(k)  [y (k)  -  x (k | k-T) ]  , 
C(k | k)  =  [ I  -  K(k) ]  C(k|k-T)  , 

where  K(k)  is  the  Kalman  gain 


(5) 


K(k)  =  C(k |k-T)  [C(k |k-T) +  B(k) ]  ]  ,  (6) 

and  B(k)  is  the  sample  noise  covariance  matrix  obtained  from  R  by 
dividing  each  of  its  elements  by  N(k)  .  The  second  equations  in  (4) 
and  (5)  are  called  the  error  covariance  equations. 


The  first  two  models  differ  in  the  number  of  observations  ob¬ 
tained.  In  Model  A  we  assume  that  the  number  of  observations  in  each 
survey  is  a  known  constant  n^(k)  =  ot'ier  words, 

P[N(k) =  nd]  =  1  . 

Model  R  assumes  that  a  constant  scalar  designated  sample  size  n^  is 

to  be  used  for  all  survey  times,  that  the  obtained  sample  size  is  a 
random  variable  with  expectation  0n^  ,  where  0<  9<  1  is  known,  and 

that  the  inter-survey  intervals  T  are  nonvarying. 
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The  t'.iird  model  (Model  C)  involves  a  continuous-time  process-- 
a  Brownian  motion  process.  Here  we  assume  that  the  sample  sizes  ob¬ 
tained  are  i.i.d.  random  variables  and  we  allow  the  intersurvey  periods 
to  vary  as  a  function  of  the  sample  sizes  obtained.  The  remaining  as¬ 
sumptions  and  the  analysis  parallel  those  of  the  first  two  models. 


tu  ail  three  models  we  are  interested  in  determining  r.  (k)  , 


or  n  , 

U 

the  designated  t 

? ample  size  to 

be  used  at  time  k  .  We  - 

Vt  rr- 

mine*  the 

n,(k)  that  minimizes  a  cost 
d 

per  unit  time  function  ; 

.  We 

define  J 

by 

J  =  E{(l/T)<c0+Clnd(k)+c2N<k»  +  L  ]  , 


where  c  is  the  fired  start-up  cost  of  ordering  a  survey, 

y/ 

unit  cost  per  designated  item-interview  (for  m  items),  and 
additional  unit  cost  per  obtained  item-interview.  Letting  a 
loss  weighting  coefficient  set  by  the  decision  maker.  L 

avg 


1 

is 


is  a 


is  the¬ 
be  a 


defined 


L 

avg 


T-l 

a  trace  [(1/T)  l  C(k+jjk)] 
j=0  ~ 

a^  trace  fC(k|k)  +  ((T-1)/2)Q]  . 


(8) 


I'he  1.  is  an  imputed  loss  attributed  to  the  use  of  survey  estimates 

of  the  state  of  the  process  and  is  proportional  to  the  sum  of  the  esti¬ 
mation  error  variances.  By  (6),  C(kjk)  is  dependent  on  the  realiza¬ 
tion  of  N(k)  through  K(k)  . 


J.  Model  A 

We  specify  for  Model  A  that  there  is  a  sequence  of  m-item  surveys 
of  one  homogeneous  population  with  obtained  sample  sizes  N(k)  =  n^(k)  . 
The  filter  gain  K(k)  of  (6)  may  be  written  as 

K(k)  =  C(kfk-T)  [C(k  | k-T)  +  R/n^k)]”1  . 
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I .et  as  now  introduce  the  concept  of  equivaLcyiU  ;:utn i-Lc  size  and  show  how 
the  error  covariance  equations  in  the  optimal  filter  theorem  may  be  use¬ 
fully  analyzed  using  this  concept.  Let  the  matrix  N  (k+jjk)  be  de¬ 
fined  by 


N.(k+j|k)  =  R1/2  C_1(k+j|k)  R1/2  ,  j  =0 , l , . . . , T- 1  , 


0 

,1/2  .1/2 


(9) 


where 

K(k)  =  R 


!\  .  The  Kalman  gain  may  then  be  written  as 
1/2  J  (  k.  j  k-T )  R1/2  [R1/2  K~X  (k  |  k-T)  R1/2  +  K/n,(k)  j-i  ,  (10) 


and  tl'.us  we  obtain  a  recursive  relation  in  N^/klk)  , 


NQ(kjk) 


nd(k)T 


Nq  (k-T  j  k-T) 


[  I  +  TR 


-1/2 


Q  R 


-1/2 


NQ (k-T | k-T) ] 


-1 


(11) 


If  only  one  item  is  surveyed  (m=l) ,  then  Q  and  R  (and  K)  are  sca¬ 
lars,  the  matrix  NQ(kjk-T)  is  replaced  by  the  scalar  Nq (k | k— T)  ,  and 

K(k)  =  nd(k)  /  [N0(k|k-T)  +  nd(k)] 

and 

N0(k|k)  =  nd(k)  +  NQ(k|k-T)  . 

This  development  leads  to  equivalent  sample  size  relations  to  replace 
the  error  variance  equations  in  the  optimal  filter  theorem  for  the 
scalar  case  of  Model  A: 

he  LWl'ch  cufVtiji; , 

N0(k+j|k)  =  NQ(k|k)  [1  +  jQR_1  N0(k|k)]_1  ; 

'it.  nup:'cy<<, 

N0(k|k)  =  nd(k)  +  NQ(k-Tjk-T)  [1  +  TQR-1  N()(k-T|k-T)]-1  . 

Here  M^(k+j|k)  is  the  equivalent  sample  size  remaining  at  time  k+j  , 

j  time  units  after  ordering  new  stock.  If  the  system  is  in  steady 
state,  Ny(kjk-T)  or  Ng(k+Tjk)  may  be  interpreted  in  inventory  terms 

Since  nd(k)  is  the  sample  size  of  the  survey 


-fe¬ 


ns  the  "reorder  point." 


conducted  at  time  k  ,  w  interpret  N  (k|k)  as  tin*  s i z»-  ox  j  .sur- 

vev  '...at  would  be  required  to  obtain  the  same  degree  of  precision  of 

estimate  at  time  k  as  that  provided  by  an  optimally  combined  data  set 

of  size  n,(k)  +  N  (k|k-T)  . 
d  0 

Figure  1  depicts  the  pattern  of  deterministic  decay  ana  replen¬ 
ishment  of  equivalent  sample  size  for  the  scalar  ease  a.  None .  A. 

Kc turning  to  tne  case  where  Q  and  k  are  matrices,  using  (S) , 
(5),  aim  (V)  in  (7)  we  obtain 

.1  (N^fk  ;  k)  ,T)  =  (l/T)  lcQ+ c^  m  n^fk)  ]  !-  a_L  trace  |  (k  |  k)  K  +  ( (T-  1 )  ,'2  )Q  ] 

=  (1/T)  [cQ  +  cx  trace  [NQ (k | k)  -  N0(k-T|k-T)  • 

(l  +  TQK_lN0(k-Tlk-T))"1]  ) 

(12) 

+  Oj  trace  [n"1 (k | k) R +  ( (T-l) /2)Q |  . 

The  optimal  pair  (n*(k)  ,  T*)  or,  equivalently,  (N*(kjk),T*)  ,  can 

easily  be  determined,  since  for  our  survey  systems  there  always  exists 

a  T  ,  denoted  T  ,  beyond  which  it  is  not  possible  to  lengthen  the 
max 

intersurvey  interval  T  without  violating  1  <  t race [ N  (k+T | k) ]  .  Thus 

we  calculate  (n*(k)  ,  T*)  by  setting  T  at  successive  integer  values, 

l’=l,2,...,T  and  computing  the  value  of  the  sample  size  n,(k)  which 

max  d 

minimizes  (  for  each  T  .  We  then  adopt  the  T  corresponding  to  the 

minimum  of  J  over  T  . 

Lt  can  be  shown  that  when  k  is  large  the  error  covariance  ma¬ 
trices  C(k)k),  C  (k+T  |  k+T)  ,  ...  ,  will  approach  a  limit,  say  C.  .  Equiv¬ 
alently,  N^(k|k)  will  approach  a  steady-state  matrix  N  .  If  the 

system  is  nearing  steady  state,  then  C(k|k-T)  wiLl  also  approach  a 
limit.  From  the  covariance  equation  in  (4),  we  then  have 
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T--+26 


C(kjk-T)  -*•  C  +  TQ 

The  steady-state  Kalman  gain  K  will  be 


so  that 


That  is , 


K  =  (C  +  TQ)  (C  +  TQ  +  R/n  )  , 


C  =  (I  -  K)  (C  +  TQ) 


TQ  =  K(C  +  TO) 


-1 


=  (C  +  TQ)  (C  +  TQ  +  R/n  )  (C  +  TQ)  , 


or 


CQ  1  C  +  TC  -  (T/nJ  R  =  0  . 

-  -v  ~  a 


(13) 


We  solve  (13)  for  C  ,  obtaining 

C  =  - ( T / 2 ) Q  +  Q1/2  [ (T/n)Q~1^ 2  R  Q~1/2  +  (T2/4) I ] 1/2  Qi/2 


(14) 


A  numerical  example  of  Model  A 

Suppose  for  m=2  the  cost  coefficients  are 

c  =  $10^  ,  c^  =  $6'/. 50  ,  c^  =  0  ,  and  =  2*10^  , 

and  the  Q  and  R  noise  covariance  matrices  have  the  following  numer¬ 
ical  entries: 


‘.00005 

.00004  ' 

and 

R  = 

'  0.  LO 

0.05  ' 

_ . 00004 

.00007  . 

.0.05 

0.25  . 

-1/2  -1/2 

Let  1)  be  the  matrix  of  eigenvalues  of  Q  R  Q 


Q  K  Q  l^2 


‘  3140 

-1978" 

then 

D  = 

'1889 

0  * 

.-1978 

5018. 

0 

1)269. 

We  can  construct  an  orthogonal  matrix  V  whose  columns  are  normalized 
eigenvectors  of  Q-1/2  R  Q~1//2  ,  so  that  Q~1//2  R  Q_1/2  =  V  D  V'  .  We 
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For  the  numerical 


will  also  have  need  ox  the  matrix  S  =  V'  Q  V  . 
values  above. 


V  = 

".8453 

-.5344" 

and  S  = 

".00009184 

.00002619" 

..5344 

.8453. 

..00002619 

.00002816. 

We  will  find  it  convenient  to  make  the  change  of  variable 
i  /n .  ,  where  0  <  n  is  the  sample  size  at  each  survey  time.  W 

Cl  u 


the  loss  function  appropriate  to  the  end-use  allocation  formula  is 
quadratic  loss 


.ssurwc 


h 

avg 


Its i ng  (141  we  obtain 


T-l 

a  trace  [  (1/T)  (C+jQ)  j  . 

j=o  -  ~ 


'avg 


a^  trace 


-0/2  +  (zl)  +  (T2/A )  I ) 1/2 


V '  0V ) 


where  we  have  used  the  property  that  tracefABC]  =  trace[BCA]  ,  that  if 
U  is  an  orthogonal  matrix  that  diagonalizes  the  matrix  A  ,  then  'd 

also  diagonalizes  A*^2  ,  or  if  A  =  UBU'  ,  then  A^~  -  UB^^h'  ,  and 
that  tracefA]  =  trace[U'AU]  =  tracelB]  ;  see,  for  example.  Bellman 
(1970). 


We  may  now  write  our  cost  function  J  for  the  two-item  repeated 
survey  problem  with  scalar  sample  size  n^  and  sampling  interval  T  as 

•J  =  (1/T)  [Cy+Cjinn^]  +  a^  trace  | -Q/2  +  ( (T/n^)  I)  +  (T2/**)  I ) 1  ^2  S] 

21/2  2  1/9 

=  c()/T  +  ct  m/z  +  aT  I  s]  i  ( z  d n  +  T  (z  d22  +  T  /4)  ~ 

~  (qll  +  q22)/21  * 

I t  can  be  shown  that  I  is  convex  In  n^  given  T  and  convex  in  T 
given  n^  .  In  practice,  (n*,T*)  will  be  found  in  a  region  in  whicli 
J  is  convex  in  (n  ,T)  .  Fixing  T  successively  at  T=1 ,2, . . . ,10 


10 


year*.,  we  .-»>>  Lve  •  ».)/  uz  -•  0  lor  z  ,  and  hence  i  ar  .<  .  -  i  / /.  ,  !>y  a 

U 

numerical  search  procedure.  We  then  compute  .1  l'or  each  T  using  the 

minimizing  n,  and  pick  the  T  and  n,  combination  with  lowest  aver- 
d  d 

age  annual  cost  .1  .  For  the  given  cost  coefficients,  numerical  results 
are  displayed  in  Figure  2.  The  T  and  combination  for  which  Liu 

average  cost  is  minimized  7=3  years  and  a  -  >,So.  u.-.i . .  .  ...... 

a  =  $875,473  per  year.  As  Figure  2  makes  clear,  tin-  average  c  <•.  i  curve 
is  rather  fiat  in  the  neighborhood  of  the  optimal  T*  and  typ icaliy  t:,«. 
survey  administrator  will  not  incur  major  additional  costs  by  choosing 
T*  -f  1  or  T*  -  1  instead  of  T*  .  It  may  also  be  seen  that  sampling 
too  frequently  is  relatively  more  costly  than  sampling  too  seldom,  as¬ 
suming  the  correctness  of  the  underlying  random  walk  process  model  for 
the  socioeconomic  variables.  An  administrator  who  is  concerned  that 
the  underlying  process  parameters  may  take  unexpected  jumps  or  exhibit 
turning  points,  which  are  not  modeled  by  the  simple  time-invariant  ran¬ 
dom  walk  models,  would  presumably  opt  for  sampling  more  frequently  than 
the  optimal  interval  found  by  this  method. 

4 .  Model  JB 

In  this  section  we  present  a  steady-state  analysis  of  a  vector 
random  yield  model  in  which  a  fixed  scalar  designated  sample  size  n^ 

is  used  for  all  survey  times,  and  T  remains  constant.  In  the  use  of  a 
fixed  order  size,  the  development  resembles  the  treatment  of  random 
supply  in  the  inventory  literature,  but  in  our  present  analysis  we  have 
level -dependent  deterministic  decay  rather  than  the  independent  stochas¬ 
tic  demands  of  Karlin  (1958).  The  stochastic  nature  of  N(k)  enters 
the  steady-state  analysis  in  the  case  of  a  fixed  n^  only  through  the 
expected  obtained  sample  size, 

E[N(k)  ]  =  0  n_,  . 

a 
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Figure  2. — Average  cost  versus  sampling,  interval:  Mode 


jii  a  random  yield  model,  an  order  inn  rule  ilnit  r.  quirts  u.->  to 
order  a  !  ix.d  amount  at  every  survey  time  is  an  "open  loop"  approach  in 
the  si  use  that  not  all  of  the  information  available  or  to  heroine  avail¬ 
able  from  the  evolving  history  of  the  survey  system  is  used  in  choosin;; 
the  sampling  interval  T  and  the  designated  sample  size  .  In  the 

ape;  fal  .  usi  i  rented  hero,  we  do  not  accommodate  tin-  genera  1  e  i.«  o.  del 


1  ll.il  I. 

■  lei. 

l  i  i  n.u; 

7  may  vary 

scocna^t i 

,  1  .  y  .  . 

-S'.  ..  . 

.  o. 

li.t.  *  i. ,  i  .i.- 

:  u..»: '  i  oa 

I'ft  Dull'S. 

.  u 

‘a  i  1 

•  1  C  ^  T  O  ,  i  * 

O  w 

.  <  -  )  i  - 

t  r.il-e  ,  l  (k 

K  W  (l.t 

'  J  /  ,  . 

\  .  i 

>  » '  t  ( r  ft  '  '  1 1 

U  1  2  it 

1  +  a,  trace 

1 F|t(k ;k) , 

r  ((V-i 

)/2  el  . 

(  •  *> 

Not  •  ;  h  4i  i  is  k  )  is  now  Si  jchast  ica  1 1  v  dependent  an  N(k.) 

sii.ee  we  as.  ume  that  the  stochastic  system  is  in  steady  si. .to  ,  vt 
I  i  ini'  ,  d  -  (,(k.k>  ;  t  1 1-  !'  k  - 1  i  ,  as  the  conditional  expei  ,  i-n  error 

i  :  a  r  i  am  e 

C  -  K|C(kjk)  i  n<1,  T,  Q,  K )  , 

where  is  a  known  constant  and  n^  =  n^(k)  for  all  k  .  Since  .) 

and  K  lire  positive  definite,  so  are  C  and  C  .  For  scalar  sample 

> i ze  n  ,  the  C  wo  seek  is  found  by  methods  similar  to  those  employed 
d 

in  Section  J  for  Model  A,  which  is  a  fixed-yield  vector  model  with 
..i  ular  siJuip  1  e  size.  The  Kalman  relations  must  now  be  expressed  in  the 
r:,\  i  ii  iutions  among  expected  error  eovar  ianees ,  expected  Kalman 
am,  and  expected  obtained  sample  size.  We  now  have 

(  =  |  I  -  K|  (C  +  TQ)  , 

'I  i  ;  e 

K  -  (C  +  U}>  |t  +  T()  +  h|K/N(k)  |  |_1  . 
this  leads  to  the  matrix  equation 

t  I)”1  t  +  1  C  -  T  K  F. j  1/N (k)  |  =  0  .  (16) 

1 1,  is  equation  is  of  the  form  of  F.quation  (13)  except  that  It  involves 
expectations  rather  than  deterministic  quantities. 
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.  -4gh 


Solving  (lb)  for  (1  wo  obtain,  strictly  paralleling  (14), 

C  =  —  (T / 2 )  Q  +  Q1/2  [T  E[l/N(k)]  Q_1/2  R  Q~1/2  +  (T2/4) I } 1/2  Q1/2  , 

so  that,  since  E(C +  ( (T-l) /2)Q]  =  C  +  ((T-1)/2)Q  ,  we  may  substitute 
the  value  of  C  found  above  into  the  cost  function  J  given  by  (15). 

The  solution  procedures  tor  Model  A  may  then  be  applied  to  1'inu  (n*,T*;  . 

) .  c*  i.  v.' 

Model  C  is  an  extension  of  random-yield  Model  B.  We  consider 
here  a  scalar  process  which  evolves  in  continuous  time  t  with  survey 
measurements  taken  at  time  points  t^  ,  k=l,2 .  Here  we  allow  the 

intersurvey  periods  T^  ,  k=l,...  (T^  =  tfc  j)  ,  to  vary  as  a  func¬ 

tion  of  the  obtained  sample  size  N(t^_^)  .  We  assume  that  the  scalar 
process  is  represented  by  the  linear  stochastic  differential  equation 

dx(t)  =  x(t )dt  +  dB(t)  ,  t(.  <  t  ,  (17) 


where  x(t)  is  a  unidimensional  Brownian  motion  process  with 


E((d6(t>)  ]  -  Q  dt  . 
The  survey  equation  (2)  now  becomes 


N(tk) 


y(tR)  =  x(Lk)  +  (1/N(tk))  l  u.(tk)  , 


(18) 


k=l ,  2 , . . .  and  VVVl’ 

where  y(tk)  and  u^(tk)  are  defined  as  before.  As  set  forth  in 

■lazwinski  (1970,  Theorem  7.1),  the  optimal  filter  for  the  system  (17) 
and  (18)  becomes  the  continuous-discrete  filter: 


[’•:thir.r>n  r.urvriir. , 

dx(t | t)/dt  »  x(t | t)  , 
dC(t|t)/dt  =  2C(t|t)  +  Q  ,  t.  <  t <  t 


(19) 


k+1  * 
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(20) 


^(tk'tk)  =  ^^k^^  +  K(tk)  ly(tk}  "  X(tk'tk)1  ’ 


C(tkitk)  =  U  '  K(tk)]  C(tk,tk)  * 


where  the  Kalman  gain  K(tk)  is  given  by 


K ( t  L )  =  C(L,  t.  )  fC(tk|tk)  +  K/K(iL)j 


(21) 


and  t  represents  the  time  instant  just  before  tlie  survey  conducted 

(and  instantaneously  processed)  at  t^  .  As  shown  in  Jazwinski  (1970), 
if  we  integrate  the  process  equation  in  (19)  over  intervals  [ t  , t^+ ^  I 
we  may  write 

tk+l 


(tk+i)  =  +  /  « 


=  *(tk)  +  w(k+l,k)  , 

where 

tk+l 

w(k+l,k)  =  /  dB(T)  . 


By  Jazwinski  (1970,  Theorem  4.1),  (w(k+l,k) )  is  a  zero-mean,  white 

(.'aussian  sequence  with 

E[(w(k+l,k))2]  =  (tR+1-  tk)Q  , 

and  the  cont i nuous-discrete  filter  for  Model  C  may  he  imbedded  in  a  dis¬ 
crete  filter  paralleling  Models  A  and  B.  The  integrated  form  of  (19)  is 


x(t|t)  =  x(tk|tR)  , 


C(t|t)  =  C(tk|tk)  +  (t-tk)Q  ,  tk<t<t 


(22) 


k+1  ' 


We  now  analyze  Model  C  with  response  rate  <j>(k)  •  Suppose  that 
the  equivalent  sample  size  on  hand  at  time  t^  is  N^(t^)  ,  The 
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.  ,  ,  ,  i  ,  ....  ,  .  .  • «  :  r.,,;ir,:n..  .it  l  i :  ^  *  l  ,  ^  ( l  t  '>  ,  .  ■"> 

1  ' 

VVV  11  1  (l-lo),',K'1  V«.  1  ■  (2,) 

Will'll  )  j  us  l  falls  to  a  prose  lotted  value  n^  (tl.e  reorder 

point),  I  n  ,  a  new  survey  of  designated  sample  size  n  ^  is  ordered 

•  ...  ■  ‘  .  .  v  .i. . :  .  r,,;  . ...  . .  .  i .  ..  • 

.  .  .  .  .  i  i  ,  ..  .  .  We  ass.i.  .e  ...  \  i  .  .  ,  - 

i 1 . .  :  .  • '  .  i  .  ,  .  ^  ,  ,  .  i . . .  J  ‘ :  i  ■ ; .  ee  ...  (t.  i  y  -  ’  .  t  -  .  J  .  .  *  ■  .ie  — 

u  :  !•  i  ; 

a  ii.se..:  .  i ;  !  i  -  a  i  \  ( .  i  ,  \  i  .  .  in  . < i  wi ,  i .:  '....,i  , 

•  v.i ■  l..ivi  i  renewal  proeess  with  renewals  oerurring  at  tiie  ;  i .a. •  points 

wi,ei  .  it  t  )  nas  just  decaved  to  the  prose  1  ee  t  ed  reorder  mint  i.  , 

i  ll  r 

"  r 


We  now  denote  bv  ,'  ( t ,  )  the 

k 

dinted  it  L,  .  We  assume  that  the 

k 

*  1 1  it  .i  known  distribution  and  that 


it)d 


i'  r« 


V'eV  ‘ 


f|V'k  v 


'  n 1 V  i  ■ 


response  rate  lor  the  sarvev  eon- 
;>(tk)  are  i.i.d.  random  variables 
0  •  v(t^)  I  .  We  tiuis  ol'lain 

n  +  : (t  )n  , 
r  k  d 

n  +  n 
r  d 


In  i.iveiitor.  term.,  we  thus  have  a  problem  of  -  toihasti.  renlen- 
. si. ,.ni,  with  designated  ordet  size  i  i  j  ,  delivered  order  size  N  ( t  ^ )  , 

,a,l  u  -da  l  ed  equivalent  sample  -  i/.e  N  ( t  ^  t  ^  )  .  d  i  veil  the  realization 

a  .ft  t  i  ,  we  then  have  level -dependent  det  erm  i  n  i  st  i  e  dee.iv  from 
0  k  l< 

( t  , t  )  down  to  the  reorder  point  n  .  I  he  time  required  for  this 
( )  K  k  r 

. . .  ■  i  .is  to  take  plane,  T  |  ,  is  found  in  terms  of  n^  ,  n^  ,  and  ‘  ( t  ^ ) 


lb 


k+i 


(k/Q)  i  1/n  -  (r>r+  'i  • 


w-  -» ) 


We  may  also  write 


k+] 


=  QK  nr  t„r  +  <Ktk)nd]  /  >(lk)nd  - 


(25) 


A  tvpii-a!  pa i  tern  of  love 1-dependent  determinist  to  doc-ay  and  eteei.  i,; 
replenishment  is  depicted  in  figure  3. 

Tin-  average  ordering  cost:,  ii  (-He,  ) )  ,  over  an  interval  oi 

avg  a 

ionjitn  l'k+1  -  tR+1  -  tR  ,  is 


,:;,VSt'!,{,k))  =  Tk+1  k0  +  (ci  +  c2<,(tk))nd] 


-] 


Cih) 


=  QK  nr  (nr  +  <Ktk)nd)  [cQ  +  (C]  +  r^d^  )»d  ]  /  <P(tk)»d 


Since  C(tit,.  )  =  K  /  N  (t|t.  )  ,  we  may  write  the  average  quadratic  loss 

K  U  K 

I.  COd,  ))  over  the  interval  T,  ,  as 
;ivj*  k  k+1 


‘k+1 


'•,v8«<'k»  -  <VW  /  C(-I.k)d-  . 


whicii  yields 


'nvB(+(,k» 


( a  ^  R / 2 )  |l/nr  +  1 / ( n  r  +  C (tk)nd)  | 


(27) 


sing,  (8)  and  (15),  we  obtain  the  cost  function  J  as 


•I  =  i:.  | c:  (4>d,  ))  +  i.  (<j>(t.  ))  ] 

4>  avg  k  avg  k 


(28) 


It  can  he  verified  by  taking  second  derivatives  that  J  is  convex  in 

n  ,  I  ‘  n  ,  for  fixed  n,  ,  0  <  n,  ,  and  in  n,  ,  0  <  n  ,  ,  for 
r  r  d  d  d  d 

lixed  n  ,  I  n^  .  Tf  we  compute  the  determinant  of  the  Hessian  of 

i  +  I,  ,  we  find  that  for  some  (but  not  all)  combinations  of  values 
avg  avg 

of  ,  n^  ,  and  the  cost  coefficients  the  determinant  is  negative, 

ilierefore,  .1  is  not  convex  in  genera].  Nonetheless,  in  practice  there 
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is  ail  aj.l  i,.:a  , 
rimvi'x  ami  i  t 

t  i  ons  l.i  /  in 


r  i  i.  (.1,0  interior  01  I'L'.-iud  in  won:.  .  i. 

mav  he  round  as  an  iterative  solution  to  the  pair  oi  e<jua- 

-  0  ,  '.).!/ •hi  -  0  . 

d 


Case  ol  Uniform  y 


u  L..O  ease 


ut  i  i  »  v.) 


u;i  , 


{  )  /  . 


•  i  i  n  .  ■+  u  ,  > 


.  /  .  I  -  A  =  (b-a)  Log (ti/a ) 


(J/n  )(i.>-a)  fog  I  (nr  1  f-'-d >  /  <nr  r  aiO 


i  L  licit 


,;i,;avK(:)|  ’  Lliav^',i 


av; 


-  i 


t  >  K  le,  n'A/n  t  e  ,  n*  A  *  e.,n“  +  on  +  on  :<  +  e.ai  aid 

1  l)  r  d  I  r  dr  Or  1  r  d  2  r  d 


(20) 


f  (a  |  K/2  )  [  I /nr  +  (  I /iitl)  (l>-a)  1  Jog|  (rir  +  lin  ^ )  /  (n  r  +  an^)  j]  . 

Now  I  mav  lie  minimized  over  (n  ,  ri  , )  by  standard  iterative  methods 

r  d 

elosely  paralleling  those1  used  for  (s,q)  inventory  problems.  A  numer¬ 
ical  example  of  Model  C  with  a  uniform  <p  is  portrayed  in  figure  A.  As 
is  siiown  in  Smith  (1980),  a  Markov  chain  set-up  may  also  lie  used  t  o  lind 
an  opt  ima 1  rn le. 
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